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We discuss the rotation of the linear polarization plane and the production of circular polariza- 
tion generated by a cosmological pseudoscalar field. We compute analytically and numerically the 
propagation of the Stokes parameters from the last scattering surface for an oscillating and a mono- 
tonic decreasing pseudoscalar field. For the models studied in this paper, we show the comparison 
between the widely used approximation in which the linear polarization rotation angle is constant 
in time and the exact result. 



PACS numbers: 95.35.+d, 14.80.Mz 

I. INTRODUCTION 

In 1977 R. Peccei and H. Quinn |l| suggested a so- 
lution to the strong CP-problem of QCD introducing a 
new symmetry breaking at a given energy scale f a . The 
boson associated with this broken global symmetry was 
called axion. All the physical properties of this pseu- 
doscalar field strongly depend on the energy scale f a at 
which the new symmetry is broken: the particle mass and 
the coupling constants with other particles are inversely 
proportional to f a . Pseudo-Goldstone bosons also arise 
in many particle physics scenarios 0j. 

Axions and, in general, other pseudoscalar particles are 
among the most favoured particle physics candidates for 
the cold dark matter (CDM) @,H M- Thc Y interact with 
photons according to the lagrangian: 



C 



(1) 



where g<f, is the coupling constant, F^ v is the electromag- 
netic tensor and = \^ vpa F pr j its dual. Many con- 
straints on axion derive from this interaction with pho- 
tons: laboratory experiments (photon-axion conversion 
experiments) and astrophysical arguments (stellar evolu- 
tion of red giants) constrain g^ to be small. One of the 
most stringent experimental bound (g$ < 8.8 x 10~ 20 
eV" 1 for m a < 0.02 eV) is obtained by the CAST exper- 
iment @ constraining the axion-photon conversion for so- 
lar axions. This limit supersedes the one obtained from 
the duration of the helium burning time in horizontal- 
branch stars in globular clusters: g^ < 10 -19 eV -1 0, 01- 
In this paper we wish to study in detail the coupling of 
such a pseudoscalar field with photons. The interaction 
in Eq. (JTJ) modifies the polarization of an electromagnetic 
wave propagating along intervening magnetic fields, or 
through a slowly varying background field <fi @ . Here we 
are interested in the second case, which does not require 
the presence of a magnetic field (note that in the first 
case the polarization is also modified in absence of axions, 
an effect known as Faraday rotation). Wc consider the 



time dependent pseudoscalar condensate as dark matter 
or part of it and study the impact of its time derivative on 
the polarization of the photons. As a consequence of its 
coupling with a pseudoscalar field, the plane of linear po- 
larization of light is rotated (cosmological birefringence) 

BO. 

In the case of Cosmic Microwave Background (CMB) 
photons, we pay attention to the rotation along the path 
between the last scattering surface (LSS) and the ob- 
server, modifying the polarization pattern generated by 
Thomson scattering at LSS This rotation induced 

by the pseudoscalar interaction modifies the gradient and 
curl of the polarization pattern (E and B following [l2l]). 
creating B modes from E modes. The parity violat- 
ing nature of the interaction generates non-zero parity- 
odd correlators (T B and E B) which would be otherwise 
vanishing for the standard Gaussian cosmological case 
\\A\ . In particular the T B power spectrum may be 
very useful to constrain the coupling constant g^ between 
photons and pseudoscalars, since it is larger than the 
auto and cross power spectra in polarization; in general, 
these non-standard correlators are already constrained 
by present data sets [H, [H, E3] • 

We study two representative examples for the dynam- 
ics of a pseudo-Goldstone field behaving as dark mat- 
ter (see [18[ for a pseudoscalar field model of dark en- 
ergy) : the oscillating and a monotonic decreasing behav- 
ior. In the latest case we study analytically the problem, 
whereas in the former numerically and analytically. The 
case of a field growing linearly in time has been studied 
in We compare the polarization power spectra ob- 
tained describing the rotation of linear polarization with 
a time dependent angle with the ones obtained consider- 
ing a constant rotation angle. 

Our paper is organized as follows. We write the rel- 
evant equations for the electromagnetic gauge potential 
coupled to a pseudoscalar field in Section II. We review 
the Stokes parameters for a monochromatic electromag- 
netic plane wave and the Boltzmann equation for CMB 
photons coupled to pseudoscalars in Section III. In Sec- 
tion IV we write the Stokes parameters in terms of the 
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left and right polarizations gauge potential and solve the 
differential equations for the latter for oscillating behav- 
ior of the pseudoscalar field. In a similar way Section V is 
dedicated to the monotonic behavior of the pseudoscalar 
field. In Section VI we test the constant rotation angle 
approximation. We conclude in Section VII. We work in 
units where the speed of light is equal to one (c = 1). 



II. ELECTRODYNAMICS COUPLED TO A 
PSEUDOSCALAR FIELD 



where k is the Fourier conjugate of z. 
These equations can be decoupled introducing 
A±{k,rj) — A x (k,ij) ± iA y (k,rf), left and right compo- 
nents of the electro-magnetic vector potential: 



A'i(k, V )+ [k 2 ±g^'k] A±(k,r l ) = 0. 



(14) 



III. STANDARD REVIEW OF STOKES 
PARAMETERS AND BOLTZMANN EQUATION 

A. Stokes Parameters 



The lagrangian density L for the photons and the pseu- 
doscalar field (f> is |2(| (following the notation of pl|): 



Iv^Wt-VM-^F^F™ . (2) 



The Euler-Lagrange equations resulting from this la- 
grangian are: 



V M F^ = 



(3) 

(4) 
(5) 



Using the definition of the electromagnetic tensor F^ v 
VM" - VM^ 1 Eq. gj) becomes: 



UA V -V V {V^-R^vA, 



'(V, 



y% pa F pa . (6) 



The complete antisymmetric tensor contain the determi- 
nant of the metric q an d [• • •] guarantees anti-symmetry 
in the four indexes [22j ]: 



-g[a/3-fS] , 



(7) 

(8) 



For a spatially flat Fricdmann-Robcrtson- Walker uni- 
verse the metric is: 



ds 2 



-dt 2 + a 2 (t)dx 2 = a 2 {r]) [-d-q 2 + dx 2 



(9) 



where t is the cosmic time, rj is conformal time and x 
denotes the space coordinates. We consider a plane wave 
propagating along h in Coulomb gauge (V • A = 0). If h 
is aligned with the z axis and neglecting the spatial vari- 
ation of the pseudoscalar field </> = 4>{rf), the two relevant 
components of Eq. ([6]) are: 



A"(r,,z) 



d 2 A x (r,,z) 

dz 2 
d 2 A v ( V ,z) 

dz 2 



= -.900' 



,9 A y (r),z) 
dz 

,dA x (r),z) 



dz 



(10) 

(ii) 



Defining Fourier transform as A X y(k,rf) = 
(2tt)^ 1 J e lkz A x y (r], z)dz the previous equations be- 
come: 

A'^k^) + k 2 A x (k, v ) + g^'ikAyik^) = 0, (12) 
A''(k, rj) + k 2 A y (k, rf) - g^'ikA x {k, rf) = . (13) 



The complex electric field vector for a plane wave prop- 
agating along z direction at a point (x, y) in some trans- 
verse plane z = zq is: 



E 



(E X (t),Ey(t)) 

e x e x (t)e^ 



e y £ y (t)e 



itp y (t) 



-ikt 



(15) 



where the physical quantity is the real part of E. For a 
spatially flat Fricdmann-Robcrtson- Walker metric the re- 
lation between the electromagnetic tensor and the phys- 
ical fields is: 



F„ 



—E x —E y —E z 

B Z —By 

B x 
-B x 




(16) 



In general we consider quasi-monochromatic waves: the 
amplitudes (e x (t) and £ y (t)) and the phases (<p x {t) and 
<p y {t)) are slowly varying functions of time respect to the 
inverse frequency of the wave. 

The Stokes parameters /, Q, U andV are defined as: 



I 

Q 

u 



-((E:(t)E x (t)) + (E;(t)E y (t))) 



((E*(t)E x (t)) (E* y (t)E y (t))) 
((E*(t)E y (t)) + (E* y (t)E x {t))) 

(S x £y COS (</? x - tPy)} , 



V ee ((E*(t)E y (t)) - (E* y (t)E x (t))) 
= -4 (e x s v sin (tp x - <p y )) . 



(17) 
(18) 

(19) 
(20) 



where (• • ■) denote the ensemble average, the average over 
all possible realizations of a given quasi-monochromatic 
wave. For a pure monochromatic wave ensemble averages 
can be omitted and the wave is completely polarized: 



I 2 -Q 2 -U 2 -V 2 = {). 



(21) 



The parameter / gives the total intensity of the radiation, 
Q and U describe linear polarization and V circular po- 
larization. Linear polarization can also be characterized 
through a vector of modulus: 



Pr. 



VQ 2 + u 2 , 



(22) 
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and an angle 9, defined as: 

. 1 



U 

— arctan — 



(23) 



It is important to underline that / and V are physical 
observables, since they are independent on the particular 
orientation of the reference frame in the plane perpen- 
dicular to the direction of propagation n, while Q and 
U depend on the orientation of this basis [23[ . After a 
rotation of the reference frame of an angle 9 (R(9)) they 
transform according to: 



Q 

u 



Rte) 



Qcos(26») + C/sin(26»), 
-Qsfn(20) + J7cos(26>) . 



(24) 



Also linear polarization, like total intensity and circular 
polarization, can be described through quantities inde- 
pendent on the orientation of the reference frame in the 
plane perpendicular to the direction of propagation of 
the wave. In the context of CMB anisotropies, the lin- 
ear polarization vector field is usually described in terms 
of a gradient-like component (E mode) and of a curl-like 
component (B mode). 

In a similar way it is possible to describe the electric 
vector field in the x — y plane through a superposition of 
left and right circular polarized waves defining: 



ie. 



and 



e_ = 



V2 



(25) 



V2 

In this new basis: 

I = ±{(E* + (t)E + (t)) + (EZ(t)E-(t))) , (26) 
Q = ±{(E* + (t)E-(t)) + (EZ(t)E + (t))) 



(i- 

2 



(e+e„ cos (tp+ -^-)) , 



U = --((El(t)E_(t))-(E*_(t)E + (t))) 
2 

= — (e+£-sin((p + -</?_)) , 
a 

V = ±((El(t)E + (t))-(E*_(t)E_(t))) . 



(27) 

(28) 
(29) 



The relation between the vector potential and the elec- 
tric field for a wave propagating in a charge-free region 
is: 



dA 

"dt 



(30) 



According to definition given in the previous section the 
Stokes Parameters in terms of the vector potential are: 



1 



' = ^((a';a' + ) + (a'za'_)) , 
q = ^{(a';a'_) + (a'ia' + )) 



(31) 



u 



V 



-± {(a';a'_) (a'ia + )) 

1^((A'*A'_)), 

i ((A'*A' + ) - (A'*A'_)) . 



(32) 

(33) 
(34) 



As we shall see in more detail in the following section, 
the coupling to a cosmological pscudoscalar field induce 
a physical time-dependent rotation of the plane of linear 
polarization along the line of sight, described by: 

Q' = 26'{rj)U and U' = -26'(r))Q , (35) 

whose solution is: 

Q = QiCos26 + Uism26, 
U = -Q l sm29 + U i cos29. 



(36) 



where Q t , C/j are the Stokes parameters at initial time 
which would be otherwise unchanged in absence of the 
interaction with the pscudoscalar field. 



B. Boltzmann Equation and cosmological 
birefringence 

In the Boltzmann equations for linear polarization of 
the radiation density contrast averaged over momenta 
contains a mixing term: 

29' = . g 0' , (37) 

due to the pseudoscalar interaction [Tlj : the Boltzmann 
equation for spin-2 functions Q ±iU is: 



A 'a±iu( k , V) + iknA Q±lU (k, rf) 



-n e a T a(r]) 



AQ±iu(k,v)+J2\[Y 



t2 KT4 m) (M) 



Ti2e'{rj)& Q±iU {k,ri). (38) 

where [i is the cosine of the angle between the CMB pho- 
ton direction and the Fourier wave vector, n e is the num- 
ber density of free electrons, ot is the Thomson cross 
section, S Y™ are spherical harmonics with spin- weight 
s, and Sp \k, rf) is the source term for generating linear 
polarization reported in [24| (m = 0, ±1, ±2 corresponds 
respectively to scalar, vector, and tenor perturbations): 



(fc, r,)+12V6A^ 2 (k, 7 ? )+12v / 6A^(fc, rj) . 



(39) 

Ay" ; and A^"; are the Fourier transforms of the coeffi- 
cients of the following series: 

A T (x,n,?/) 

= ^(-z)V^(2/ + l)A^ ) (x, 77 )rr(fi) , (40) 

/ m 

AQ±ic/(x,n,?7) 

= £H)VM2Z + i)Ai m )(x, n) ±2 vr (fi) , (4i) 



4 



Note that Eq. (|38|) corrects some typos in Eq. (1) of Ref. ^ B = _1 2 ^ +iU _ $ 2 A Q ^ (43) 

[l8| |. 2 \ / 

The quantity Aq±w is related to the rotation invariant 
polarization fields A^ and through the spin raising 
(d ) and lowering (d ) operators: 

Y , - s Following the line of sight strategy for scalar pertur- 

A^; = — — (d 2 Aq + hj +d 2 Aq-hjJ , (42) bations we obtain, in agreement with Ref. (l8| : 

I 



A T (fc,?7o) = / dr)g(r))S T (k,r))jt(kr) -kr)), (44) 

J Vice 

A E (k, m ) = r d V g( V )S { °\k 1 r 1 ) J f T]0 - fc g cos [2%)] , (45) 
Jv™ (kr]o - kr)) 

A B (k, m ) = r drjg^S^ik, ^ 710 - ^ sin [20(7?)] . (46) 
Jfhaa (kijo - kr)) 

where g(rf) is the visibility function, ST(k,rj) is the sorce term for temperature anisotropics, and ji is the spherical 
Bessel function. The polarization Ci auto- and cross-spectra are given by: 

Cf E = (4^) 2 ^2), / k ' dk [A*(Mo)f , (47) 

Cf B = ( 4 -) 2 L^ +2 l I k ' dk [As(fc,%)] 2 , (48) 



16 (e-2)\ 

9 (£ + 2)1 
16 (£-2)! 



C f B = ( 4 ^) 2 TF7T 1 % I k 2 dkA E (k,r) )A B (k,r) ), 



Cj E = (4tt) 2 J A / fc 2 dfcA T (Mo)A B (Mo), (50) 



Cf 3 = (4^) 2 J jL^LJl I k 2 dkA T (k, Vo )A B (k, Vo ). (51) 

I 



9(£- 


1-2)! 


16 (£ 


-2)! 


9(f- 


1-2)! 


16(1 


-2)! 



In the approximation in which 6 = 9, with constant 
in time, Eqs. f4"B")) and (|4"6"|) simplify since cos[2(9] , sin[20] 
can be extracted from the integral along the line of sight 
and: 

A°i s = A E (6 = 0) eos(20), (52) 
A°^ s = A E {6 = 0) sin(20) , (53) 

and the power spectra are given by [Til . fl5| : 



n EE,obs 


= Cf E cos 2 (26), 


(54) 


n BB,obs 

u e 


= Gf E sin 2 (26), 


(55) 


n EB,obs 
U i 


= l -C EE sin(A6) , 


(56) 


n TE,obs 
U £ 


= Cj E cos(29), 


(57) 


^TB,obs 


= Cf E sm(29). 


(58) 



The expression for 9 to insert in Eqs. (52-58) is: 

= f ttM - MVrec)} ■ (59) 



Several limits on the constant rotation angle 9 have been 
already obtained using current observation of CMBP (see 
Tab. H]). 



Data set 


9 (2a) [deg] 


WMAP3 and Boomerang (B03) [15] 


-13.7 < 9 < 1.9 


WMAP3 [16] 


-8.5 < 9 < 3.5 


WMAP5 [17] 


-5.9 < < 2.4 


QUaD [25J 


-1.2 < 6» < 3.9 



TABLE I: Constraints on linear polarization rotation 9 in the 
constant angle approximation. 

This time independent rotation angle approximation 
is an operative approximation which allows to write Eqs. 
(f53")) . is clearly inconsistent since for 9 = const the term 
proportional to 9' in the Boltzmann equation (|38[) van- 
ishes and therefore there is no rotation of the linear po- 
larization plane. See Figs. [5] and [7] for a comparison 
of this approximation with a full Boltzmann description 
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of the birefringence effect for a dynamical pscudoscalar 
field. 



IV. COSINE-TYPE POTENTIAL 

In this section we assume that dark matter is given by 
massive axions, <fi is governed by the potential [3[: 



V{4>) = m- 



f 1 

J a 

N 2 



cos 



4>N 



(60) 



where N is the color anomaly of the Peccei-Quinn sym- 
metry. Here we are interested in the regime where the 
axion field oscillates near the minimum of the potential 
(for simplicity we shall consider N = 1 in the following) : 
4>l fa 1 and the potential can be approximated with 
V{4>) — m 2 cj> 2 /2. In this case (f>(t) satisfies the equation: 



3H<f> + rn 2 



0. 



(61) 
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FIG. 1: Evolution ofH./7i. ICC in function of conformal time for 
m = 10" 28 eV (dashed line), m = 5 x 10" 27 eV (dotted line) 
and for a matter dominated universe (continuous line), from 
recombination (r) re c) to 3.5r7 rcc . Present time corresponds to 
?7o = TfcccV 1 + Zrcc — 33.1877 rcc . 



When ra > 3H the scalar field begins to oscillate, and 
the solution in a matter dominated universe (a/ a = 2/3i) 
is (Hi: 



<f>(t) = t~ 1/2 [ Cl J 1/2 (mt) + c 2 J_ 1/2 (mi)] 



mt>l (p0 . r \ 

~ — sm [rat) , 
rat v ' ' 



(62) 



where the time-independent coefficients of the Bessel 
functions c\ , c 2 depend on the initial conditions. 

The averaged energy and pressure densities associated 
with the field are: 



1 n— 2 mt3>l u 



-p- cb 2 1 

Pa. = m 

2 2 



2 
n_ 

2t 2 
2 

2t 2 



[ — + 

rat 



2~I 2 m ^ >1 W. x 0| _ 



1 

mi 



,(63) 
(64) 



where ~ denotes the average over an oscillation period of 
the axion condensate. Note that we are implicitly assum- 
ing that the pseudoscalar field is homogeneous. In the 
context of axion physics, this means that in our observ- 
able universe we have just one value for the misalignment 
angle, which means that the PQ symmetry has occurred 
before or during inflation. 

We fix the constant 4>q comparing with the energy 
density in a matter dominated universe: 



Pm 



3# 2 M 2 



67rf 2 



pi 



3tt 



Afpi 



3irrat 



sin(?7it) , 



(65) 



(66) 



where M pl ~ 1.22 x 10 la GeV is the Planck mass. 



Using the relation between cosmic and conformal time 
in a universe of matter: 



3 \Vo 



we find the following approximation for (f>(rf) 
4>(rj) — \l — Pl . 1 sin 



and 



rar]o 
3" Mpi 

7T 77 



Vo V 



(67) 



(68) 



cos 



V0 / 7? 



— ^ sm 



varj 



>)() 

Vo ( V_ 

Vo 



.(69) 



If m is not too small the value of Ti. = a' /a obtained with 
the scalar field density in the Friedmann equation coin- 
cides with that of a matter dominated universe H. = 2/rj 
once the average through oscillations is performed [27j 
(see Fig. [TJ . The derivative can be replaced in Eq. (fl4|) 
for the evolution of the gauge potential: 

A'i(k,r 1 ) + k 2 [1 ± A( mg ^m,k, m )} A ± (k, v ) = 0, (70) 

defined the function: 

'3 <? Af pl 



A(j7;00,m,fc,77o) 



ir kr\ 

H . 

t sir 

mr\ 6 



cos 



Vo ( V_ 



.Vo ( V_ 

no 



(71) 



This term, induced by axion-photon coupling, oscillates 
with frequency proportional to the mass of the axion and 
its amplitude decreases with time. 
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In the next two subsections we study analytically and 
numerically Eq. (fT0|) for different values of the param- 
eters m and we exclude the region where the mass 
of the pseudoscalar field is so small that the field starts 
to oscillate after equivalence (m < 3H cq ), and the re- 
gion corresponding to a PQ symmetry broken at energies 
higher than Planck scale (f a > M p \): see Fig. [2] 



A. Adiabatic solution 



Adiabatic solutions of Eq. (|70|) are: 



A, 



1 



a ±i / uo s dri 



(72) 



1 ± 



g<p<t>'(y) _ 



fcyl ± A(?y) and s 



where Lu s {r]) = k\ 
±. 

The second derivative respect to conformal time is: 



A„ — A s 



4w 2 



The adiabatic solution ([72")) is a good approximation for 



3w .; 



the vector potential when the terms t— 
small compared to w 2 : 



and 



3A' 



<M 16/c 2 (l±A) 3 
^ _ ±2(1 ± A)A" - A / 2 
2w? ~ 8fc 2 (l± A) 3 

If both condition are satisfied and A <C 1: 



< 1 



r are 

(74) 
(75) 



A± 



exp 



^2fc(l±A/4) 
1 

^2fc(i±7r 50 ^fy 



»7± 2 



Adry 

exp [±z (fc?? ± 2tt 5 0)] . (76) 



In the adiabatic regime the coupling between photons 
and axions produces a frequency independent shift be- 
tween the two polarized waves, which corresponds to a 
rotation of the plane of linear polarization: 

^adiabatic = y [<t>(Vo) ~ HVrcc)} ■ (77) 

This result agrees with the one obtained in Ref . [1[ , which 
therefore holds in the adiabatic regime. More important 
than this, ^adiabatic = i-c. Eq. (|77[) agree with the 
rotation angle which is approximated by Eq. (|59|) in the 
Boltzmann section III.B. This agreement is not a coin- 
cidence and shows the usefulness of studying the gauge 
potential as done in this section: the estimate based on 
the adiabatic approximation of the rotation angle due to 
cosmological birefringence can be also obtained by study- 
ing the gauge potential A s . 



Typically <p(rj ICC ) 3> 4'{ r lo)\ from last scattering to now 



p ~ m 



3 AfpiWfa) 



3 M p i / rio 



8ir m V 27T rnrjo \ r\ 

An estimate of the angle ^adiabatic is: 

Tm p1 



'adiabatic — 9<fi^ 



(1 + z rcc ) 3/2 - 1 



(78) 



(79) 



8tt mrjo 

Note the dependence of ^adiabatic on the coupling con- 
stant and on the mass of the pseudoscalar field: for fixed 
g<j, the effect is larger for smaller masses. 

The amplitude of the electromagnetic field changes ac- 
cording to: 

2 



E 



A' 



LPs 

2t? ' 



(80) 



so the degree of circular polarization evolves according 
MM: 



U c = 



A' 



A' 



A' 



At 



yO + A - VT^A _ A 
v/T+A + v/l"! ~ 2 



(81) 



B. CMBP constraints on the (m,g^) plane 

In a flat universe dominated by dust (w = 0) plus 
a component with w = — 1 (cosmological constant) the 
evolution of the scale factor in terms of cosmic time is 
M: 



a(t) = 



n 



MAT 



i - n 



MAT 



sinh 



■yl — ^MATHot 



(82) 

where J^mat is the density parameter for matter nowa- 
days. The Hubble parameter is: 



H = H y/1 - MAT coth [ - Vl - SlMArHot ) . (83) 
The pseudoscalar field evolves according to: 



[sinh (|V1 - S W#o*)] 



x sin 



mt\ 1 - (1 - Omat) 



3ffo 
2m 



(84) 



The energy density is 



1 2T 2 



mt»l 



2[sinh(§Vl-fiMAT^oO 



oc a" 3 . (85) 
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Assuming that the axion-like particles contribute to the 
cold dark matter density p < j >i o = ^mat Pck, o (where 
PGR, o is the critical density) we can estimate fa: 



3(1 — SImat) H Q M pl 
7r 2m 



(86) 



Therefore the evolution of the pseudoscalar field as a 
function of cosmic time is: 



3^MAT H M p i 

n 2ma 3 / 2 (i) 



x sin 



mt\ 



(1 - Omat) 



V 2m 



.(87) 



Note how this equations reduces to Eq. (|66|) in a matter 
dominated universe: SImat = 1, H /(2a 3 ^ 2 ) = 1/ (3t). 
The linear polarization plane, from last scattering sur- 
face, rotates according to: 

m = ^ im - 0(w)] . (88) 

The Boltzmann equation contains the derivative of 
the rotation angle respect to of conformal time (cfr. 
Eq. ([38]) ). so we need the relation between cosmic and 
conformal time. For a particular model with I^mat = 0.3 
it is possible to fit numerically the relation between cos- 
mic and conformal time from last scattering to nowadays: 



Vo 
3.45 



3.09 



(89) 



Replacing this expression in Eq. lj5T)) we obtain the evo- 
lution of the pseudoscalar field as a function of conformal 
time 4* = 4>(rj). 

The linear polarization angle is not constant in time, 
but it oscillates with varying amplitude. If the field rep- 
resents a fraction J^mat of the universe energy density, 
then the amplitude of these oscillations is: 



^(^MATjTO,^) 



1 J 3Qmat g(j>MpiH ( 1 
4' 



1 



3/2 



3/2 



351mat g4>Mp\Ho ^ 3j/2 



(90) 



Fixed f^MATj it is possible to constraint a certain re- 
gion of the (to, <70)-plane requiring ^(^mati m , 9<f>) to 
be smaller of a certain angle, typically of the order of few 
degrees (see Tab. UJ. The excluded region considering 
current limits on CMB birefringence is shown in Fig. [2] 
Fixed a particular value for the pseudoscalar field mass 
and for its coupling with photons we can also estimate 
how the polarization angular power spectra are modi- 
fied by a rotation of the linear polarization plane. We 
modified the source term for linear polarization in the 
public Boltzmann code CAMB [3l[ following Eqs. (gSJ) 
and (|46p . The linear polarization rotation angle is given 
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FIG. 2: Plane (log 10 m [eV] , log 10 [eV" 1 ]): region ex- 
cluded by CAST [f| (blue with vertical lines), region where 
|#a(Omat = 0.3, m, g$)\ > 10 deg (red region with horizontal 
lines), (m, g^) values expected in main QCD axion models 
(red with dots), region where the mass of the pseudoscalar 
field is too small in order to explain dark matter (m < 3-ff oq ) 
(yellow with horizontal lines), and region where PQ symme- 
try is broken at energies higher than Planck scale (f a > M p i) 
(yellow with vertical lines). 



by Eq. (|88|) and the evolution of the pseudoscalar field 
by Eq. (|87[) . The new power spectra are compared with 
the standard unrotated ones in Fig. [3] fixed to = 10~ 22 
eV and g<p = lO-^eV" 1 . 

In Section IVII we compare the power spectra modi- 
fied version of CAMB obtained starting by Eqs. (|45I46|) 
which takes into account the time dependence of the 
pseudoscalar field in the integral along the line of sight 
with the approximated spectra obtained following Eqs. 
(55-59). 



C. Comments for axion cosmology 

For axions the coupling constant with photons and 
the energy scale f a at which the new symmetry is broken 
are related 0]: 

W = g^ with o.i<e<i, (91) 

where the value for £ depends on the particular model 
considered for the axion. By using this relation a limit 
on the coupling constant is turned into a limit on the 
energy of symmetry breaking. 
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FIG. 3: EE (a), BB (b), TE (c), TB (d), and E B (e) angular power spectra for m = 10" 22 eV and gj, = 10" 20 eV _1 (black 
solid line), the black dotted line is the standard case in which there is no coupling between photons and pseudoscalars (9 — 0). 
For the B B power spectrum (b) we plot for comparison also the polarization signal induced by gravitational lensing (black 
dotted line), and primordial BB signal if r = 0.1 (blue dot-dashed line). The cosmological parameters of the flat AC DM model 
used here are Q, b h 2 = 0.022, O c h 2 = 0.123, r = 0.09, n s = 1, A s = 2.3 x 10~ 9 , H = 100 /ikms" 1 Mpc" 1 = 72kms _1 Mpc -1 . 
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The critical density associated with the misalignment 
production of axions strongly depends on the initial 
misalignment angle associated with the axion field 0^ 
through the following relation [H, Q| : 



0.23 x 10 



±0.6 



fa 



10 12 GcV 



1.175 



e 2 F(e s 



(92) 

where h encodes the actual value of the Hubble param- 
eter (Hq = 100/ikms -1 Mpc -1 ) and F(Oj) accounts for 
anharmonic effects if 0j ^> 1. The demand fi m j s < SIdm 
provides an upper bound on /q' 175 © 2 (assuming F (0j) ~ 

i) Hiiniag: 



f a e}- 7 < 2 x io 



11H-12 



GeV. 



(93) 



This condition becomes also an upper bound for f a 
under the assumption that inflation occurred before the 
breaking of PQ-symmetry (f a < /inf) 01- m this sce- 
nario different regions have different values for Oj, so 
averaging over all observable universe the value of 0j in 
equation can replaced by its rms value (tt/v3) and the 
limit f a < lO 11 ^ 12 GcV is obtained. As can be seen from 
Fig. 2, CAST disfavors values of g 4 ~ lO" 1 *-- 12 GcV -1 
with a mass up to 0.02cV. Note however that our cal- 
culation cannot be applied directly to this case since we 
assume (f> homogeneous in our universe, whereas it is not 
if the PQ symmetry breaking occurs after inflation: al- 
though taking into account space inhomogeneities were a 
second order effect in cosmological perturbation theory, 
cosmological birefringence might be larger than the one 
computed in this paper. 

Our calculations apply without modifications to the 
case in which inflation occurs after PQ-symmetry break- 
ing: the initial misalignment angle G^ is homogeneous 
throughout our universe and can be much smaller than 
n/y/3. Such possibility allow the scale of PQ-symmetry 
breaking f a to be much higher than lO 11 ^ 12 GeV and is 
motivated by anthropic considerations [HI, HH, US H3- 
These smaller values of can be constrained by present 
data in CMB polarization in a much better way than 
CAST, in particular for small masses. 



V. EXPONENTIAL POTENTIAL 

We consider in this section a pseudoscalar field with 
an exponential potential: 




FIG. 4: For A = 4.5 Dashed line: J7dm +fibaryons, dotted line: 
Orad, thin continuous line: fi</>, thick continuous line: w^, in 
terms of the natural logarithm of the scale factor (from log a ~ 
—40 to nowadays logao = 0). Here ODM.o + ^baryons.o = 0.852 
and JVo = 0.148. 



about a vanishing value (see Eq. (f69|) ), we study here a 
case where the behaviour is monotonous. 

It is known [H[ that exponential potential with A 2 > 
3(1 + wf) leads to a component which tracks the domi- 
nant background fluid with equation of state = w^p^. 
In order to satisfy the nucleosynthesis bound we choose 
A = 4.5. During the matter dominated era the scalar 
field behaves as: 



+ Vo exp (-Ak</>) = /pmat = f PMA ^° (95) 



P<t> = -y - Vo cxp (-\nq 



(96) 



where PMAT = PDM + Pbaryons + P4>- 

For A = 4.5 the contribution of the pseudoscalar field 
to universe energy density is shown in Fig. [4] The value 
of f2^ changes with time, but it is almost constant (f2^ ~ 



< !> .0 



0.148) from recombination (loga r , 



-7) to 



nowadays. 

The derivative of the pseudoscalar field respect to con- 
formal time is proportional to a -1 / 2 and the evolution 
of the scale factor in the matter dominated phase is 
a(v) = (n/vo) 2 so: 



n 770 

V J PMAT,0 — 
V 



(97) 



V ((/>) = Vo exp (— Ak<5 



(94) 



with K 2 = 8ttG . Theoretical motivations to this kind of 
potential are certainly weaker than the ones for the po- 
tential presented in Eq. (|6"0|) . However it is interesting 
to show how the kinematics of the pseudoscalar field is 
important for the resulting spectra of CMB anisotropics 
in polarization. Whereas the time derivative of the pseu- 
doscalar field in the previous case contains oscillations 



Substituting this relation in Eq. (|T4|) we obtain the fol- 
lowing expression for the evolution of the electromagnetic 
potential: 



k 2 ± g+y/fpuATfi — k )A±=0 



(98) 



This is a particular differential equation, called Coulomb 
wave equation; defining q± = T5*x//Pmat,o»7o/2 = TO. 



10 



and x = kr\ it becomes: 



dx 2 



x 



-4-1 



(99) 



The solution of this particular equation can be written 
in terms of regular (Fq (Qjx)) and irregular (Go(q,x)) 
Coulomb wave functions [33, [|(| : 

A+ = f+F (q + ,x) + g+G Q (q + ,x) 
= f+Fo(-q,x)+g+G (-q,x), 

A- = f-F (q-,x)+g-G {q-,x) 
= f-F (q,x)+g-G (q,x), 



where /+, /-,.<?+,.<?- €E C; in a compact notation: 
A±{q,x) = f±F (Tq,x) + g±G (Tq,x) ■ 



(100) 



The Stokes parameters contain the derivative respect to 
conformal time 77, so we evaluate: 



A' ± (q,x) 



dF (Tq,x) dG (Tq,x) 

j± 7; \-g± 



dx 



d. 



(101) 



The solution given in Eq. (jlOOp verifies the Wronskian 
condition (A±A± — A'±A± = i) if the following relation 
holds: 



f±g± - f±g± = T 



9(/ifl±) = ^. (102) 



In the general case, when the coupling does not van- 
ishes (t/0 7^ 0), we expand the solution (|100ll for large 
value of x neglecting terms proportional to 0(x~ 2 ) (see 
Appendix) : 



A±(q,: 



f± |^ cos (x ± a (q, x)) + (l T ^) sin (x ± a (q, x)) 

2 

+g± [ 1 T cos (x±a (q, x)) - — sin (x ± a (q, x)) 



where a {q, x) = q In 2x — arg T(l + iq). The derivative respect to conformal time is: 

,,2 



A' ± (q,x) 



1 ± — cos (x ± a (q, x)) 
2xJ 



2x 



sin (x ± a (q, x)) 



-g± 

k 
' 2 



— — cos 

2x 



(x ± a (q, x)) — ^li -^-^j sin (x ± a (q, x)) 
(/± 



„i(a:±a(«,s)) f x ± JL , • 1. 

2x 2x 

„2 



-i(x±a(q,a;)) 



i± 



2,r 



*^ ' (/± ~ i9±) 



r 



(103) 



(104) 



In general both forward moving waves (A± cx e lkv ) and 
backward moving waves (A± cx e lfcl? ) must be taken into 
account for propagation of light in a medium. Chosen a 
particular value for the constants f± and g± that verifies 
the Wronskian relation (|102p the evolution of polarization 
is fixed. 

If we assume, according with [4l], |42j , that the photon 
pscudoscalar conversion is a small effect due to low en- 
ergy of CMB photons, the production of backward mov- 
ing waves can be neglected (see [43| for the use of this 
approximation). The Eq. (|104[) setting f± = —ig± be- 
comes: 



A'± (q, x) 



2x 2x ,' 



(105) 



and in terms of the value at recombination time 

q fl 1 



A'± (q,x) ~ A' ± (<7,x rec ) 



1± 



.9 









*^rec / 



cxp{-i [ 

& ^rec 

± Aa]}(^06) 



where we have introduced 

Aa = a{q,x) - a(q,x iec ) = gln(r//?7 rec ) 
= -ln(o/ore C ) ■ 



(107) 



We observe that also in this exact case the plane of linear 
polarization is rotated of an angle Aa independent on k 
whose dependence on the difference between the present 
value of 4> and the corresponding one at recombination is 
the same of the adiabatic approximation and of Eq. (|59[) . 
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FIG. 5: EE (a), B B (b), TE (c), T B (d) and EB (e) angular power spectra for gj, = 10" 28 eV" 1 (black solid line); the black 
dotted line is the standard case in which there is no coupling (0 — 0). For the B B power spectrum (b) we plot for comparison 
also the polarization signal induced by gravitational lensing (black dotted line), and primordial BB signal if r = 0.1 (blue 
dot-dashed line). The cosmological parameters of the flat CDM model used here are flj, = 0.0462, Q. c = 0.9538 (fi<p ~ 0.148), 
t = 0.09, n s = 1, A a = 2.3 x 10" 9 , H = 72kms~ 1 Mpc" 1 . 
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Current measures and constraints on the polarization 
pattern of CMB anisotropics produce an upper limit on 
the linear polarization rotation angle of the order of few 
degrees (see Tab. fl}. We now use these constraints and 
our analytic expression: 



The additional term in the Boltzmann equations which 
rotates the linear polarization plane is (see Eq. 



Ti26'(r])A Q±iU (k,r]), 
which clearly vanishes for 8 = const. 



(Ill) 



|0| = ^ln(l + z rec ) 



- JY^^,o|^|M p iln(l + 2 rec ), (108) 

to obtain an upper bound for q, which can be turned into 
a upper bound on g^; if \6\ < 6 deg, then: 



M^lCr^eV- 1 



(109) 



where we have assumed: Sl^.o — 0.148 and z roc ~ 1100. 

The angle of linear polarization 8(rj) appearing in 
Eqs. (|4"S]) and l[4"6]) can be replaced with: 



l%)l 



1 



3 / i] 
i „— ^0,0 MM pl ln 



(110) 



and the polarization power spectra arc evaluated using 
the expression given in Section IIII B[ see angular power 
spectra of Fig. [5] 

In Section I VII we compare the power spectra modi- 
fied version of CAMB obtained starting by Eqs. (I45l46j) 
which takes into account the time dependence of the 
pseudoscalar field in the integral along the line of sight 
with the approximated spectra obtained following Eqs. 
(55-60). 



VI. COMPARISON WITH CONSTANT 
ROTATION ANGLE APPROXIMATION 

In this section we compare the angular power spec- 
tra obtained modifying the public code CAMB [3l[ con- 
sidering the correct dynamic of the pseudoscalar field 
(8 = 8(rj)) as described in Section Till B[ with the ones 
obtained in the constant rotation angle approximation 
(8 = const) for the two different potential considered in 
the previous sections: see Figs. [6] and [Jj 

In Section [III Bl we have already shown how the power 
spectra in the constant rotation angle approximation 
[Eqs. ([54]) - ([55]) ] can be obtained from the general expres- 
sions [Eqs. |47|) - ([5T|) ]. 

Power spectrum modifications obtained starting di- 
rectly form the Boltzmann equations and taking into ac- 
count the temporal evolution of the pseudoscalar field 
are usually smaller than effects predicted considering a 
constant rotation angle equal to the total rotation an- 
gle from last scattering to nowadays. If the cosmological 
pseudoscalar field evolves quickly then the constant rota- 
tion angle approximation clearly leads to an overestimate 
of the effects. 

It is important to stress that the constant rotation 
angle approximation is just an operative approximation. 



VII. CONCLUSIONS 

We have studied the impact of a pseudoscalar field act- 
ing as dark matter on CMBP. We have shown that such 
pseudoscalar interaction with photons rotates the plane 
of linear polarization and generates circular polarization. 
In absence of measures for the V mode of CMBP, the 
existing upper limits on an isotropic T B and E B corre- 
lations can constrain the coupling constants of photons 
with the pseudoscalar field. 

We have examined two representative examples for the 
dynamics of a pscudo-Goldstone field behaving as dark 
matter: the oscillating and the monotonic decreasing be- 
havior. In the monotonic decreasing behavior, by ne- 
glecting backward moving waves, we have shown how 
present CMB observations can constrain the coupling 
constant to small values as O (10~ 30 ) eV. For the more 
physically motivated axion case which leads to an oscil- 
lating behaviour, we have shown how constraints from 
CMB cosmological birefringence can become important 
for small masses for the axion. 

We have also shown how the use of integral solution 
of the Boltzmann function may improve the estimate ob- 
tained by multiplying the CMB power spectrum by the 
suitable trigonometric functions of the rotation angle as 
in Eqs. ^MB- 
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Coulomb wave equation 



The Coulomb wave equation is [3 

d 2 w 
dx 2 



1 2q L(L+1) 



w = 0, 



(112) 



with x > 0, — oo < q < oo, L a non negative integer. 
Here, in order to solve Eq. ([99]) . we are particular inter- 
ested to the particular case when L=0. 

The solution can be written in terms of regular 
(-Fi(<7,x)) and irregular (GL(q,x)) Coulomb wave func- 
tion: 



w = ciF L (q,x) + c 2 G L (q,x) . 



(113) 
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FIG. 6: EE (a), B B (b), TE (c), T B (d) and E B (e) angular power spectra for m = 10~ 22 eV and 94> = 10~ 20 eV _1 
(black solid line) and approximating the rotation angle with the constant value 8 rcc (red dashed line). The cosmological 
parameters of the flat AC DM model used here are Q b h 2 = 0.022, Q c h 2 = 0.123, r = 0.09, n s = 1, A s = 2.3 x 10~ 9 , 
H = 100/ikms" 1 Mpc" 1 = 72kms" 1 Mpc" 1 . 



The Coulomb functions can be expanded for large val- similarly for the first derivative respect to x 
ues of x [39| : 

F = gcos8 + fsin6, (114) F„ = g* cos (9 + /* sin (9 , (116) 

G = /cos(9- ff sin<9, (115) G = /* cos(9 - g* sin(9 , (117) 
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FIG. 7: EE (a), B B (b), TE (c), T B (d) and EB (e) angular power spectra for = lCT 28 eV -1 (black solid line) and 
approximating the rotation angle with the constant value 9 lac (red dashed line); the black dotted line is the standard case in 
which there is no coupling. The cosmological parameters of the flat CDM model used here are fij, = 0.0462, tt c = 0.9538 
(fi^ ~ 0.148), t = 0.09, n s = 1, A, = 2.3 x 10^ 9 , H = 72kms- 1 Mpc -1 . 



with 6 = x — qhi2x + argT(l + iq) and: .„ /•* * 

oo oc 



5fe 

fc=0 A;=0 
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where: 



fS 
= 1 - 

a k 



/o = 1 , fk+l = dkfk - bk9k 

9o = , 9k+i = a k g k + b k f k 

fk+l 



0> fk+i = a kf k ~ b k g* k 



q 

x 



9k+i = a k9*k + b k f k 

,2 



X 

9k+i 



2(k + l)x ' 
Restricting to the first order 

/ = 1 



r 



(2k + l)q , q 2 -k(k + l) 

h - 



2(k + l)x 





VO 


f 1 


2x ' 


yX 2 


l_ 


VO 


t 1 


2x 


yX 2 




VO 


t 1 


2x 


X 2 



(118) 
(119) 

(120) 

(121) 

(122) 



(123) 
(124) 
(125) 



(126) 



Summarizing the asymptotic expansion of Fl (q, x) and 
for large values of x is: 



F (q,x) ~ — cos( 



1 + ^ ,Sm( 



G (q,x) 
and for the first derivative 



1 + — ) cos < 

2a; / 



sin ( 



2a; 



sin t 

2a; 



G^z) 



2a; 



cos 61 - 1 



2.r 



sin . 



(127) 
(128) 
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